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^ . ABSTRACT 

< 

OO I NGC 4013 is a nearby Sb edge-on galaxy known for its "prodigious" Hi warp and its "giant" 

tidal stream. Previous work on this unusual object shows that it cannot be fitted satisfactorily 
^^ , by a canonical thin-|-thick disk structure. We have produced a new decomposition of NGC 4013, 

r) I considering three stellar flattened components (thin-|-thick disk plus an extra and more extended 

component) and one gaseous disk. All four components are considered to be gravitationally 

O ,, coupled and isothermal. To do so, we have used the 3.6/im images from the Spitzer Survey of 

JL I Stellar Structure in Galaxies (S'^G). 

^ ■ We find evidence for NGC 4013 indeed having a thin and a thick disk and an extra flattened 

C/2 , component. This smooth and extended component (scaleheight zec ^ 3 kpc) could be interpreted 



as a thick disk or as a squashed ellipsoidal halo and contains '^ 20% of the total mass of all three 
stellar components. We argue it is unlikely to be related to the ongoing merger or due to the 
off- plane stars from a warp in the other two disk components. Instead, we favor a scenario in 
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which the thick disk and the extended component were formed in a two-stage process, in which 
an initiahy thick disk has been dynamicahy heated by a merger soon enough in the galaxy history 
to have a new thick disk formed within it. 

Subject headings: galaxies: individual (NGC 4013) — galaxies: photometry — galaxies: spiral 
— galaxies: structure 



1. Introduction 

Thick disks, first detected by Burstein (1979) and Tsikoudi (1979), are seen in edge-on galaxies as ex- 
cesses of light a few thin disk scaleheights above the galaxy mid-planes. They are known to be ubiquitous 
(Yoachim & Dalcanton 2006; Comeron et al. 2011a) and their properties give important clues for under- 
standing galaxy formation and evolution (Comeron et al. 2011b and references therein; hereafter COllb). 
Recent studies (Robertson et al. 2006; Elmegreen & Elmegreen 2006; Brook et al. 2007; Richards et al. 2010) 
suggest an in situ formation mechanism for a significant fraction of the thick disk mass. This was recently 
supported further by our result that the masses of thick and thin disks are of the same order (COllb). 
Other formation mechanisms such as disk kinematical heating due to its own overdensities (Villumsen 1985; 
Hanninen & Flynn 2002; Schonrich & Binney 2009; Bournaud et al. 2009) and the accretion of stars from 
infalhng satellites (Statler 1988; Gilmore et al. 2002; Abadi et al. 2003; Navarro et al. 2004; Martin et 
al. 2004; Read et al. 2008) also contribute to the thick disk mass (COllb). 

In COllb we made thin+thick disk decompositions of 46 edge-on galaxies using images from the Spitzer 
Survey of Stellar Structure in Galaxies (S^'G; Sheth et al 2010). However, two galaxies, ESO 079-003 
and NGC 4013, could not be successfully fitted down to low surface brightness due to the presence of an 
extra light component not accounted for in our fits [affecting significantly the luminosity profiles starting at 
M3.6Ami(AB) = 23magarcsec~^]. In COllb, NGC 3628 also presents an obvious third component, but at a 
lower surface brightness [/i3.6Aim(AB) — 24.5 mag arcsec"^]. The aim of this Letter is to study the properties 
of these extra components. As ESO 079-003 is located near a bright star which makes its study at low 
surface brightness levels difficult and NGC 3628 presents a disturbed morphology due to a recent interaction 
we focused on NGC 4013. 

NGC 4013 is an Sb galaxy (Buta et al. 2007), at 18.6±2.5Mpc (NED average for 13 redshift-independent 
measurements). Its optical diameter is D25 = 294" (HyperLEDA; Paturel et al. 2003). It is often described 
using superlative adjectives: it has a "prodigious" Hi warp (Bottema et al. 1987; Bottema 1995; 1996) and 
a "giant" stellar tidal stream with an age of a few Gyr (Martinez-Delgado et al. 2009; hereafter MD09). The 
Hi warp starts just at the optical edge of the galaxy (Bottema 1995) and it is one of the largest warps ever 
observed. It may have been triggered by the minor merger event which caused the tidal stream (MD09). In 
addition, NGC 4013 has a boxy bulge (Jarvis 1986) which could be tracing an edge-on bar (Athanassoula & 
Misiriotis 2002; Martinez- Valpuesta et al. 2006) or be the consequence of a merging event (Binney & Petrou 
1985). 

The S'*G 3.6/im-band image used for the study presented in this Letter is shown in Fig. [TJ 
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Fig. 1. — 3.6/ini-band S^'G image of NGC 4013 with two different luminosity stretches. The vertical lines 
indicate the limits of the bins for which luminosity profiles have been produced, from left to right, at 
galactocentric distances of — 0.8r25 < R < —0.5r25, — 0.5r25 < R < — 0.2r25, — 0.2r25 < R < 0.2 r25, 
0.2 r25 < R < 0.5 r25 and 0.5 r25 < R < 0.8 r25. In order to avoid the influence of the bulge we have not 
produced fits for the central bin. 



2. Fitting procedure 

In COllb, the observed luminosity profiles were fitted with synthetic profiles resulting from coupling 
two stellar plus one gaseous isothermal disks (using the equations in Narayan & Jog 2002). The solutions 
are not analytic, so the fit was done by comparing the luminosity profile to a grid of pre-computed models 
with different thick-to-thin central mass density ratios (pTo/pto) a-nd thick-to-thin velocity dispersions in 
the vertical direction ratios (crT/ct) where T stands for the thick disk and t stands for the thin disk. Both 
observed luminosity and pre-computed profiles were scaled to have a mid-plane luminosity equal to unity 
and to have p{z = 2QQ)/ p{z = 0) = 0.1. In addition the synthetic profiles were convolved with a Gaussian 
kernel with a full width at half maximum (FWHM) equal to 2.2" in order to account for the point spread 
function of the 3.6/im-band image (S^'G "super-PSF"; Sheth ct al. 2010). The mass-to-light ratios (T) for 
both thin and thick disk stars had to be assumed in order to make the computation of the grid of models, 
and reasonable limiting cases with Tx/Tt = 1.2 and Tx/Tt = 2.4 were studied (see COllb). The gaseous 
disk was assumed to have a column mass density equal to a 20% of that of the thin disk. The case without 
gas was also studied yielding similar results to the with-gas case, but the fit quality was slightly worse. 

The fits were made ignoring the mid-plane pixel and over a range of magnitudes, A?n, which was defined 
as the range for which either the square root of the y^ of the fit was smaller than yx^ < 0.1 magarcsec^^ or 
down to a limiting magnitude p,\ ( AB) = 26 mag arcsec"^ . The discussion on how we chose this goodness-of-fit 
parameter and this limiting magnitude appear in COllb. 

In COllb, for the NGC 4013 luminosity profiles at 0.2 r25 < \R\ < 0.5 r25 the fits could only be 
performed down to /i(AB) ~ 23 magarcsec"^, at which point our criterion for the quality of the fit, ■\/x^ > 
0.1 magarcsec"^, was not met. The reason for that was the presence of a third component which made our 
two stellar disk decomposition insufficient at describing NGC 4013. This component was called "halo" by 
van der Kruit & Searle (1982) and by MD09 who reported that it is box-shaped (although the boxiness is 
probably caused by the location of the tidal streams and not by the intrinsic shape of the component). 

We adapted the code used in COllb to include a third component in equilibrium with the two stellar 
disks and the gaseous disk. This is a natural approach, because the luminosity profile decays exponentially 
with height once in the range of heights for which the more extended component dominates, and because the 
solution of the equations for a set of coupled flattened components can be approximated as exponential at 
high distances above the mid-plane. However, the fact of adding a third flattened component does not imply 
that it is a disk; this is because, without kinematical information, a squashed ellipsoidal halo nature for this 
component cannot be discarded. Arguments in favor of this component being a disk are (i) its relatively 
high surface brightness, and (ii) its luminosity profile parallel to the mid-plane direction (not shown here) 
which is that of a typical edge-on disk, with a shallow slope for low R and an exponential slope at larger 
R. Arguments in favor of the extended component being an squashed elliptical halo are (i) its isophotes 
are not disky (MD09) and may be very close to elliptical if the light of the tidal streams could be removed, 
and (ii) its ellipticity (e = 0.63 when measured in an ellipse fit between radius 120", which is the truncation 
radius of the thin disk, and 170", where the extended component starts to be highly affected by noise) is 
compatible with that of simulated elliptical haloes (see, e.g., Lee et al. 2005). Due to its uncertain nature 
we will hereafter term this component as "extended component" (EC). 

The grid of models we used for the fits had pto/pw, Pto/peco, crT/ct and ctt/cec as free parameters, 
where EC denotes the extended component. Our grid of models has been computed including a gaseous disk 
and using the same normalizations as in COllb. We considered the cases Tx/Tt = 1.2 and Tx/Tt = 2.4. 
In both cases we set Tx/Tec = 1-0. This is justified by the fact that once above the mid-plane dust lane 



(z > 10", a region mainly influenced by the thick disk and the EC), the colors remain roughly constant, 
as can be seen in Fig. [21 Fig. [5] provides extra evidence that the thick disk and the EC contain old star 
populations since the average color for z > 10" is g' — r'{z > 10") = 0.73 mag, which is compatible with the 
average colors predicted for SO galaxies {g' — r' = 0.68 mag; Fukugita et al. 1995). The reddening for z < 10" 
is caused by the mid-plane dust lane. Since we found no significant influence of this mid-plane dust on the 
3.6/xm-band profile for this galaxy in COllb, we considered the mid-plane dust lane to have no effect on our 
fits. 



3. Results 

The addition of a third stellar component allows us to fit the surface brightness profile down to /i(AB) = 
26.0 mag arcsec^^. The fit is significantly better than that published in COllb (Fig. |3] and Table [T]) in the 
sense that it goes far deeper. The fact that for 0.2 r25 < \R\ < 0.5 r25 Am is 3.0magarcsec~^ larger than 
what was achieved in COllb and that the outermost fitted z goes from 25" — 45" in COllb to 70" — 100" in 
the present work argues for the correctness of the fit, and thus for the physical reality of the third component. 
The numbers exclude a possible alternative, namely that the better fit is achieved merely by the addition of 
more free parameters. 

The fitted properties of the thin and the thick disk are fairly similar to those reported in COllb (Table[T]), 
and the ratio of the column mass densities Ex/St is equal, within 10%, to that obtained for the three bins 
for which we succeeded to obtain a fit in COllb. The mid-plane stellar densities are not much affected by 
the inclusion of the EC, and its main effect is to reduce the scaleheights of the thin and thick disks. The 
small effect introduced by the EC in Ex/St is due to the fact that the mid-plane density of the EC is ^ 10 
{^ 10) times smaller than that of the thick disk and ^ 20 (^ 10) times smaller than that of the thin disk 
for Tx/Tt = 1.2 (Tx/Tt = 2.4). The fact that possible undetected ECs do not have much infiuence on the 
results in COllb shows the robustness of our fitting approach in general. 

Each of the three stellar components has a significant range in vertical height for which it dominates 
the luminosity profile, making it easy to determine its properties and degeneracies unlikely. The thin disk 
luminosity dominates for z < 5", the thick disk dominates for 5" < z < 20" and the EC dominates for 
z > 20". The only problem in the fit is that the thin disk scaleheight is poorly constrained due to poor 
sampling (zt is around 1.5", which is smaller than the FWHM). 

The EC, if considered to be a disk, has an exceptionally large scaleheight, zec ~ 3 kpc, comparable only 
to that of thick disks in NGC 0678, NGC 4437 and NGC 4565 in COllb. The scaleheight is significantly 
larger for the — 0.8r25 < R < — 0.5r25 bin probably because its luminosity profile is affected by the brightest 
loop of the tidal stream (MD09). The pto/peco ratio for 0.2 r25 < \R\ < 0.5 r25 is significantly higher (a 
factor of two) than it is for 0.5 r25 < \R\ < 0.8 r25 implying that the scalelength is significantly larger than 
that of the thin and the thick disk. 

Using the weightings in Eq. 5 from COllb we find that the ratios of masses of the stellar components are 
Afx/Mt = 1.22 (Mx/A/t = 2.19) and Mt/Mec = 2.74 {Mt/Mec = 2.69) for Tx/Tt = 1.2 (Tx/Tt = 2.4). 
When taking the mass of the thick disk and the EC together (Afx + Mec) /Mt = 1.67 ((Mx + Meg) /Mt = 
3.00). With a circular velocity speed of Vc = 181.7kms"^ (HyperLEDA), NGC 4013 would not fit in 
the Mt/Mi — Vc relationship discovered by Yoachim & Dalcanton (2006) and shown in Fig. 12 of COllb 
both when considering M^ /Mt and (Mx -I- Mec) /Mt. Furthermore, if we consider the thick disk as part 
of the thin disk and the EC to be the only thick disk in NGC 4013, we get Mec/ {Mt + Mx) = 0.20 
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Fig. 2. — g' — r' color profiles made for 0.2 r25 < \R\ < 0.5 r25. The profile has been produced using SDSS 
DR7 data (Abazajian et al. 2007). Points for z < 40" are calculated every three pixels (~ 1.2") and those 
with z > 40" are calculated every nine pixels (~ 3.6"). The error bars represent the 2a statistical errors. 
The vertical lines separate the regions dominated by the luminosity of the thin disk, the thick disk and the 
EC for Tx/Tt = 1.2. 
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Fig. 3. — Example of a fit to the vertical luminosity profiles of NGC 4013 considering Tx/Tt — 1.2. Data 
points with error bars (2cr statistical errors) represent the observed luminosity profile, and the dashed green 
curve the best fit. The dotted red curves indicate the contributions of the three stellar components. The 
dash-dotted vertical lines indicate the limits of the range in vertical distance above the mid-plane used for 
the fit. The horizontal line represents the /i(AB) = 26.0 mag arcsec"^ level down to which the fits have been 
done. The inset shows the thin-f-thick disk best fit obtained for NGC 4013 in COllb. 



(Mec/ {Mt + Mr) — 0.26), which is too low to fit into the M^/Mt — Vc relationship. 

The luminosity profiles could not be satisfactorily fitted by a sum of two sech (z/zq) functions, but were 
acceptably fitted by a sum of two exponential functions. The •y/(x^) was, however, worse than in the three 
disk fit and misses a substantial amount of light corresponding to the thinner component of the three-disk 
fit (0.2 — 0.3magarcsec~^ for the inner bins). This fit, although simpler than that made with three coupled 
stellar components, should be regarded as unphysical. 



4. Discussion 

Just as causal links have been sought between the "prodigious" warp and the merging event causing the 
"giant" stream, one may also be tempted to find links between a minor merger and the EC of NGC 4013. 

If the EC is linked to the minor merger, then it can be made of tidal debris. According to the images in 
MD09, the tidal loops should not affect the luminosity profiles at i? = 0. Although producing a fit for the 
minor axis using our procedure is impossible due to the presence of the bulge, we produced the luminosity 
profile (Fig. |4]). It is clear that in the bin with — 0.2r25 < R < 0.2 r25 the EC is present and probably 
slightly brighter than in the 0.2 r25 < \R\ < 0.5 r25 as predicted in the case of a disk or an ellipsoid. The 
color analysis made by MD09 shows that the tidal stream is much redder than the EC, although the large 
error bars make a common origin still possible. Other arguments against a tidal origin for the EC are its 
smoothness, its symmetry and the uniformity of its scaleheight with varying radial distance R. 

Another possibility is that the EC has been formed by the dynamical heating of the disks by the crossing 
of the dwarf galaxy causing the tidal stream. However, as the disk self-gravity which counteracts the disk 
heating is higher at low galactocentric radii, disks formed in this way appear fiared and, as a consequence, 
the EC would form with a relatively small scaleheight at low R. Simulations (Quinn et al. 1993; Walker et 
al. 1996; Kazantzidis et al. 2008; Bournaud et al. 2009) show how this mechanism could reasonably produce 
the observed scaleheights at high galactocentric radius, but would fail to produce such an EC for the inner 
kiloparsecs. Additionally, it seems difficult to create a component containing 20% (26%) of the mass of 
the galaxy if we assume Tx/Tt = 1.2 (Tx/Tt = 2.4) in a recent event without disturbing the disks. The 
regularity and the symmetry of the galaxy argues for an old origin for the EC, thus confirming that the 
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Fig. 4. — Luminosity profile for the — 0.2r25 < R < 0.2 r25 bin (solid line) compared to that obtained from 
the 0.2 r25 < \R\ < 0.5 r25 bins (dotted line). 



assumption of equilibrium is reasonable. Finally, the stream in NGC 4013 is Monoceros-like (MD09), thus 
probably a very minor merger (^^1:100; Peiiarrubia et al. 2005) not likely to cause large disturbances. The 
possibility of the EC being created during an older merger event cannot be discarded; Purcell et al. (2010) 
have shown that a 1:10 merger 5 Gyr ago in a Milky Way- like galaxy yields results which qualitatively match 
the appearance of NGC 4013 except for the fact that their resulting EC is more boxy than the one we detect 
and we do not find the same flaring. They show that low-latitude accretion events could heat up to 1% of the 
disk mass to z > 7kpc, which is in rough agreement with what we found (1.8% of disk mass at z > 7kpc). 

We also discard the EC to be the effect of a stellar warp seen in projection. Bottema (1995) found that 
the stellar disk is truncated at the radius at which the warp starts and the warp model made by Bottema 
(1996) shows that its line of nodes is in the direction of the line of sight, making a stellar warp look as a 
thick disk when seen in projection impossible. Even if the Hi warp and the stellar warp were decoupled 
and with different lines of nodes the stellar warp would be unlikely to mimic an EC unless the stellar warp 
started at a radius smaller than the H i warp. 

Thus, we deduce from our findings and from literature data that the EC is a real feature and not some 
tidal feature or warp seen in projection. We also deduce that both the thick disk and the EC are old. 
A possible scenario for the formation of a three component system would be having a disk formed thick 
with stars formed before and during the build-up of the galaxy from small fragments (Robertson et al. 2006; 
Brook et al. 2007; Richards et al. 2010) and/or a disk heated by the internal thickening caused by kinematical 
heating due to giant clumps (Elmegreen & Elmegreen 2006). Then a merger event would further thicken 
the disk, which would become what we know as the EC. After the merger, the canonical thick disk would 
form with the same mechanisms as the EC formed prior to the merger event. Finally, the remaining gas, 
plus that coming from cold flows would settle in the mid-plane and form the thin disk. This disk formation 
mechanism is not very frequent, since only two galaxies over 30 exhibit a bright EC in COllb. The difference 
in the formation mechanisms may explain why these two galaxies fall outside the Mx/Mt — Vc relationship 
discovered by Yoachim & Dalcanton (2006). 



5. Conclusions 

NGC 4013 is not adequately described by the canonical thin-f-thick disk description (COllb). In this 
Letter, the luminosity profiles of NGC 4013 are fitted satisfactorily using the solutions of three stellar 
flattened components plus one gaseous disk in equilibrium. The newly described extended component (EC) 
has a relatively low surface brightness, but due to its vertical extent, contains a significant fraction of the 
disk mass (between 20% and 26% depending on the assumed Tt/Tx). The EC has a longer scalelength than 
the galaxy disks, is smooth and its properties do not depend strongly on varying galactocentric radii. 

The nature of the EC is unknown and could be a second thick disk or a squashed elliptical halo. The 
smoothness of the EC makes it unlikely to be related to the ongoing minor merger of NGC 4013. We also 
discard the EC to be made of off-plane stars of a warped disk, since the warp has been modeled to have 
its line of nodes in the direction of the line of sight (Bottema 1996). We favor a scenario in which the EC 
was formed in a two-stage process, in which an initially thick disk was dynamically heated by a merger soon 
enough in the galaxy history to have a new thick disk formed within it. 
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Table 1: Result of the fits. 



Fitting range 


Tx/Tt = 1.2 


Tx/Tt = 2.4 




yUi = 25.96 mag arcsec"^ 




Am = 6.0magarcsec~^ 




PTo/pto ==0.35 


PTo/pto = 0.64 




Pto/pecq = 7.00 


Pto/peco = 6.40 




CTx/o't = 2.55 


cx/c^t = 2.69 


-0.8 r25 <R< -0.5 r25 


ctt/o-ec = 0.38 


ctx/ctec = 0.38 




Sx/St = 1.21 


Sx/St = 2.23 




St/Sec = 1-54 


Sx/Sec = 1.46 




zt = 13Gpc 


zt = 120 pc 




zt = 610 pc 


zx = 580 pc 




ZEC = 3660 pc 


ZEC = 3500 pc 




yUi = 25.99 mag arcsec"^ 




Am = 7.5magarcsec~^ 




Pto/pw = 0.36 


PTo/pto = 0.60 




Pto/peco = 12.00 


Pto/peco = 12.00 




O'^jot — 2.45 


CTx/fTt = 2.69 


-0.5 r25 <R< -0.2 r25 


ctt/o-ec = 0.41 


ctx/ctec = 0.41 




Sx/St = 1.19 


Sx/St = 2.13 




St/Sec = 2.91 


Ex/Sec = 3.02 




zt = 110 pc 


Zt = 110 pc 




ZT = 510 pc 


zx = 530pc 




ZEc ==2780pc 


ZEC = 2890 pc 




p\ = 25.93 mag arcsec"^ 




Am = 7.5magarcsec~^ 




PTo/pto = 0.35 


PTo/pto = 0.64 




PTo/PECO = 11.60 


Pto/peco = 10.60 




cTT/cTt = 2.60 


CTx/cTt = 2.69 


0.2 r25 <R< 0.5 r25 


ct/ctec = 0.44 


crx/o-EC = 0.43 




Sx/St = 1.25 


Sx/St = 2.25 




St/Sec = 3.13 


Sx/Sec = 2.94 




zt = 110 pc 


Zt = 100 pc 




ZT = 520 pc 


zx = 500pc 




ZEC = 2540 pc 


ZEC = 2470 pc 




pi = 25.67 n 


lagarcsec"^ 




Am =: 5.5 n 


lagarcsec"^ 




PTo/pto = 0.39 


PTo/pto = 0.70 




PTo/pEco = 5.57 


pto/peco = 5.00 




(Tx/CTt = 2.35 


CTx/cTt = 2.50 


0.5 r25 <R< 0.8 r25 


o't/o'ec = 0.44 


CTx/o-EC = 0.44 




Sx/St = 1.19 


Sx/St = 2.19 




St/Sec = 1-58 


Sx/Sec = 1.48 




zt = 150pc 


Zt = 150 pc 




ZT = 620pc 


zx = 590pc 




ZEC = 2860 pc 


ZEC = 2700 pc 



Note. -Tx/Tt: ratio of the thick and thin disk mass-to-light ratios. [i\: limiting magnitude of the fit. Atti: 
dynamical range over which the fit has been produced. PTt)/ pto, Pto/ Peco'- niid-plane thick to thin and thick 
to EC mass density ratios, crx/crt, ctx/ctec: thick to thin and thick to EC vertical velocity dispersion ratios. 
Sx/St, Sx/Sec^ thick to thin and thick to EC stellar column mass density ratios, zt, zx, zec^ thin, thick 
and EC scaleheights. 



